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We derive a condition for the existence of completely or exponentially localised Majorana bound 
states (with the potential for non-Abelian statistics) in a generic many-body system. We discuss 
the relationship between the existence of these operators and the protection of the ground state 
degeneracy from local perturbations. We use our methods to study the exponential decay of the 
Majorana bound states in the non-interacting Kitaev chain, finding complete agreement between our 
many-body calculation and single-particle results. We then apply these results to various interacting 
systems which have previous evidence for Majorana bound states. 


In 1936 a real solution to the Dirac equation was dis¬ 
covered by Ettore Majorana [T]. This corresponds to 
a class of fermions that are their own anti-particles. 
Though physical examples of these Majorana fermions 
have eluded us to this date, the zero-energy Majorana 
bound state (MBS) has been proposed as a condensed 
matter quasiparticle analogue. This has generated a large 
amount of interest due to the suggestion it may obey non- 
Abelian statistics ElE], with implications for both funda¬ 
mental physics and quantum computation [4j. Over the 
last few years, the set of hypothetical systems supporting 
Majorana bound states has grown in leaps and bounds 
[ME]- Experiments have demonstrated the existence of 
a zero-energy state in superconducting wires piin] and 
atomic chains m, and schemes have been developed to 
implement and observe braiding experimentally nans]. 

In the last few years, interest has developed in inter¬ 
acting systems that support MBSs [HHii]. Evidence 
for topologically nontrivial excitations has been found 
in Josephson currents miiiH], entanglement spectra and 
correlation functions mini 121 ], or by continuously tun¬ 
ing to a non-trivial topological phase in a non-interacting 
limit m uni- Importantly, in some of these systems, 
strong interactions are required for the appearance of 
these effects. However, in the single-particle formalism, 
data about topological protection and the potential for 
non-trivial braiding may be read from the solutions to 
the single-particle Hamiltonian. In the many-body case, 
this data is obscured. 

In this work, we present a method for determining this 
data for a generic many-body system. We derive condi¬ 
tions for the existence of an operator that we claim cor¬ 
responds to a localised MBS. We demonstrate how the 
form of this operator may be derived when these con¬ 
ditions are satisfied, and discuss the implications of the 
existence of this operator on ground state protection. We 
then apply these methods to the previous studied models 


of [HHII]- 

For the purposes of this paper, we will concern our¬ 
selves with fermionic systems with discrete sets of sites, 
and assume all Hamiltonians conserve fermion number 
parity. Then, for a zero-energy Majorana bound state, 
we require a degeneracy between the lowest even particle 
number (I^'b)) and odd particle number (Idto)) many- 
body states of our system. To simplify the discussion, 
we will assume no additional degeneracies exist; though 
these are necessary for non-Abelian statistics, our meth¬ 
ods immediately generalise. Given such a system with 
Hamiltonian PL, suppose we wish 7 to be an operator 
corresponding to the creation or destruction of an MBS. 
We could require 7 to be Hermitian and unitary, to excite 
I^'e) to e^^\A>o), and to commute with PL These 

conditions are satisfied by the MBS in non-interacting 
systems such as m 1 ng. However, these conditions 
assume that excitations in the system behave indepen¬ 
dently of the state of the system, which is not necessarily 
the case when interactions are present. This can be seen 
in the Hubbard model in the strong-coupling limit: at 
half-filling the system is a Mott insulator, but the addi¬ 
tion or removal of a single electron turns the system into 
a conductor, with a dramatically different energy spec¬ 
trum. As such, we restrict our conditions to applying 
only on the subspace spanned by our degenerate ground 
states. This is equivalent to assuming that the gap be¬ 
tween these states and those at higher energy is much 
greater than the temperature. The above conditions then 
become 


7 '*' = 7, 

(1) 

('I'a 7^ = ^AB, 

(2) 


(3) 

{'SA\[n,^]\^B) = o, 

(4) 


where A and B are either O or E, and the sign on the 
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phase shift is dependent on whether we are going from 
|d>B) to l^'o) or vice-versa. 

These conditions can be simplified somewhat. Equa¬ 
tion 1^ is a consequence of equation and the ground 
states being degenerate. Then, equation implies that 
we can write -|- |^) with = 0, 

and expanding equation with A = B = O gives that 
1 ^) = 0 . That is, 7 is still required to strictly excite be¬ 
tween and I'f'o). This also makes physical sense, 
as otherwise 7 would not be an excitation between our 
ground states. Finally, with this new restriction equation 
[^is automatically fulfilled. We then have a set of reduced 
Majorana conditions 

y'i' = 7, ( 5 ) 

71 ' 1 'a) = (6) 

with A being the opposite parity to A. 

If 7 satisfies these, we say it is a Majorana operator. 
However, if every Majorana operator corresponded to an 
MBS, then an MBS could be found whenever there exists 
a ground state degeneracy, for the following Majorana op¬ 
erator always satisfies the reduced Majorana conditions: 

70 = cos(0) (|'I'o)('E'e| + l'E's)('I'ol) 

-I-isin((/)) (|d;'o)(4'E| - |4 'e)('I'o|) ■ (7) 

Indeed, any operator of the form -\- G, with G a self- 
adjoint operator that does not act on our ground states, 
will be a Majorana operator (and every Majorana oper¬ 
ator in our system is of this form). We call the bare 
Majorana. 

This result by itself is thoroughly uninteresting, and is 
definitely not unknown. Clearly the conditions currently 
stated are not sufficient to guarantee the system can at 
all be described as supporting MBSs with exotic prop¬ 
erties. As such, we turn to these properties for further 
restrictions. We consider the protection of the ground 
state degeneracy, and the possibility for braiding statis¬ 
tics. 

The notion of braiding of excitations is intrinsically 
tied with locality. Given two excitations, we wish to de¬ 
fine a non-trivial closed path that these excitations can 
travel in real space such that the system obtains a non¬ 
trivial geometric phase upon completion. This cannot be 
done unless the location of the particles is well-defined. 
But, if we attach the geometric phase to a creation opera¬ 
tor - c’l’ —)■ e*“c^, then the corresponding number operator 
n = c^c remains invariant. This implies that if our Ma¬ 
jorana operator contains products of number operators 
with single creation and annihilation operators, only the 
latter must be localised to allow braiding. 

Let us fix a region of our system R that we wish to 
localise our Majorana operator 7 within, and let be 
the set of sites outside R. We require that if j is a site 
in J, 7 contains no single products of Cj or cj. This is 


equivalent to requiring that 7 commutes with fij. Then 
7 will commute with every product fij = 
number operators from subsets J G jj. Any matrix ele¬ 
ment of the commutator must then be zero, and 

we calculate 

('I'ol [fij.G] = 0 

(^'olhjyl'I'^;) = {'^o\inj\'^E) 

K'l'ol^jl'f'o) — ('I'_E|hj|'I'£;)| := dj = 0. (8) 

We call dj the number difference on J, and the require¬ 
ment for dj = 0 for all J C the localisation condition 
for our system to admit a Majorana operator 7 localised 
away from (or to R). Amazingly, although we have 
currently only proven that this condition is necessary for 
a Majorana operator, it is sufficient also. In the suppli- 
mentary material, we prove this by giving a method of 
constructing the Majorana operator, including explicit 
formulae for small systems. 

Let us now consider the effect of a perturbation V upon 
the Hamiltonian. We assume that V continues to con¬ 
serve Fermion number parity, as breaking this will cause 
quasiparticle poisoning. Then, V cannot couple our two 
ground states (as they are of different parity), and so to 
first order in perturbation theory, the induced gap be¬ 
tween the ground states is 

A= |(^'£;|F|4-£;)-(d'o|I^|^'o)|. (9) 

This shows that the Majorana condition on J is equiv¬ 
alent to the requirement that our system is protected 
against small local perturbations of the form V = 6 nj. 

This is not the only type of local perturbation that 
can occur in our system. Thus, the presence of localised 
Majorana operators as defined do not entirely protect the 
system. To get more insight into this statement, we can 
insert a Majorana operator 7 into equation 

A = |(4'o|[7^7-^]|^o)|. (10) 

We see that an equivalent condition for the protection of 
the degeneracy is the existence of some (not necessarily 
localised) Majorana operator 7 that commutes with V. 
In a system where the Majorana condition holds, but the 
ground state is not protected against all local perturba¬ 
tions, an experiment to observe the effect of the MBSs 
in the system will need to be isolated from these unpro¬ 
tected perturbations. 

The above discussion implies that a non-trivial system 
must be able to localise Majorana operators away from 
all local regions (i.e. dj = 0 for all local J). This makes 
sense, as our degeneracy supports two anti-commuting 
Majorana operators, constructed from the bare Majo- 
ranas 70 and 7.n-/2 (note that it is impossible with only 
one degeneracy to find more than two mutually anti¬ 
commuting Majorana operators). When these operators 
can be spatially separated, any local perturbation that 
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FIG. 1. Measurement of the exponential decay of the Majo- 
ranas in the Kitaev chain via the ground state number dif¬ 
ference di (black line), and via the solution to the single¬ 
particle Hamiltonian (blue line). Although the single-particle 
data only exists in the non-interacting regime, the di scaling 
is measurable in any scalable many-body system. 


commutes with either one or the other cannot lift the 
degeneracy without first pushing the operators together. 

Now, let us consider a non-local region R in our sys¬ 
tem such that each subregion S C R satisfies ds = 0 , but 
da 7^ 0 . Then, we can consider i? to be a ’Majorana- 
supporting subsystem’; every Majorana we write must 
have some density in R. Note that such an R must exist, 
although it may be the entire system. When R is not the 
entire system, there is the possibility for Majorana oper¬ 
ators to either take the forms ju (if there is no number 
difference outside i?), Jr + Jr or Jr'Jr (or some combi¬ 
nation thereof) where R is the region outside R. In the 
first situation, the MBS can be considered to be split into 
parts, whilst in the second situation, the MBS will be an 
effective description of the system inside R. 

Within each region R however, it is not at all obvious 
where the MBS itself lies; we may remove all but any 
one site within the region. As an example, consider the 
Kitaev chain with Hamiltonian 

N-l N 

'^Kitaev — ^ ^ h.C.) -f /T ^ ^ Tlj . 

i=i 1=1 

( 11 ) 

in the special point where A = t and /r = 0, the following 
Majoranas are localised to every site 1 < j < N: 

f-i 

7j-0 = cos{<p)ic] + Cj) n(l - 

N 

-b isin(0)(c] - Cj) (l-2h/). 

1 = 1+1 

This is clearly a problem with our description. To rectify 
it, we consider the experimental probes that might be 


used to determine the MBS position. These usually take 
the form of some transport measurement; we connect the 
system to a lead which breaks local parity conservation, 
and attempt to locally excite between the two ground 
states. Then, we are looking for a probe density pr = 
('I'£;|pr.|'I'o) ^ 0, with pr an operator entirely confined to 
a small region r that depends on the system and probe in 
question, but with no restriction to satisfy the Majorana 
conditions. We expect to see exponentially localised 
as it is in the Kitaev chain with p^. = cl for r = {i}. 

It may well be asked here why this is not a sufficient 
condition in the first place. However, we stress that we 
can only find Majorana-like excitations when the localisa¬ 
tion condition is satisfied for all local regions. Thus, a full 
study should include demonstrating that all Majorana- 
supporting subsystems are non-local, and then investi¬ 
gating both degeneracy protection and the Majorana lo¬ 
cation. 

In a realistic system, we expect an MBS to have ex¬ 
ponentially decaying tails, and so all MBSs will exist ev¬ 
erywhere (and braiding will have an exponentially small 
chance of errors). This should cause a breaking of the lo¬ 
calisation condition in local areas, but only by an amount 
that is exponentially small in the system size, as this 
measures the minimum density of all MBSs on a given 
site. To check this, in figure we plot the deviation 
di = |('k£;|hi|'I'£;) — ('l'olhil'I'o)l for the Kitaev chain 
as we increase the system size N, in the topologically 
non-trivial phase (Hamiltonian given in equation II with 
t = 0 . 7 A, p = O.IA). Comparing this to the density of 
the single particle Majorana operators in a larger chain 
(50 sites) finds complete agreement up to a constant fac¬ 
tor. As this procedure can be repeated for a many-body 
system of adjustable size, we have a method of measuring 
the exponential decay of Majoranas in a realistic system. 
However, as noted before, this procedure by itself will 
not locate individual MBS within the system. 

We now apply our procedure to the interacting Kitaev 
chain studied in EHIIII]- Taking the formalism of [ 20 j . 
the Hamiltonian of the system is 


L-l 

H = TiKitaev + U ^ (2hj — I)(2hj_|_i — I), (12) 

i=i 


where L is the length of the chain. A ground state de¬ 
generacy is obtained when the bulk chemical potential is 
p = p* := A -\-tU-\- — A^)/ 4 )^^^, and the edge 

sites have a chemical potential half this size m- At this 
point, we can calculate the number difference and the 
probe density, with probe Pn = cJ): 


(I + a2)L ’ 

■(I-a^)"-! (l-a2)L-n 

“ [ (I + a 2 )n “ (l + a2)i-"+i 


(13) 

(14) 
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Here, a = [^*(-\/(2A/^*)2 + 1 + 1)/(2A)]^/^. We see the 
expected exponential decay of both quantities. 

We now apply our method to the double quantum dot 
of [ISl [T^ . The effective Hamiltonian of this system is 


U = -Ez + 




i=i 


+ t ^cos (^0 c\ ,^C 2 ,a + cr sin c\ ,^C 2 ,s + h.c.^ 


cr=± 


+ Ae* = cos ( — ) I ^ c\^c\^ + 2 ^ 


A At 


- sm - c „c. 


At At 


+ h.c.) 


(15) 


In the infinite U and vanishing Ez limit, there is a de¬ 
generacy when A = [T^ . The ground states are then 

given by 


l^is) = 2 - 6"*'^+/^ cos(6»/2)(c|^c^^ -h 4t4;) 

+e”*'^+^^sin(6»/2)(c4c^^-hc44i) l^)> (16) 

l^o) = ^ cos(6»/4)(cl^ - c^) 

-hsin(6»/4)(cl^-hc4) |u). 

(17) 


The description of the MBSs in this system is compli¬ 
cated by the addition of spin. If we redefine our spin 
axes by writing = cos{pj/2)c]^ + a s\n{pj/2)c]^, the 
localisation condition is satisfied with Pj = f + (~l)'^ f ■ 
That we can only do this is indicative of a lack of pro¬ 
tection in the system. Specifically, the system is not 
protected against on-site magnetic fluctuations save in 
the pj direction. The system is protected against non¬ 
magnetic fluctuations of the form fij^ + fij^. However, 
as a Zeeman potential Ez ~ ksT is required to lift the 
problematic Kramers degeneracy, a spin-orbit coupling 
of near 0 = tt is required so that this does not break 
the MBS-supporting degeneracy also. At this point, the 
probe density ('I'£;|4o-l'I'o) = independent of j and 
O’, and we can write down exact Majorana operators due 
to the small size of the system (details in the supple¬ 
mentary material). For example, a Majorana operator 
localised to site 2 with (() = (()+ = 0 is 

7 = %(;^%% + liiT%)(4t + ^2t) 

+ %(;^%% - + %)■ 

Finally, we consider the system of m- This consists 
of two spinless wires a and 6, with creation operators a\ 



Number of Sites 

FIG. 2. (top) Demonstration of exponential decay in the dou¬ 
ble wire of m via the number differences di (black, solid) and 
02 (blue, dashed). Red line is an exponential fit to 02 for ease 
of viewing, (bottom) plot of probe density as a function of 
system size, both in the middle and at the edges of the sample 
(probe defined in text). 


and b], and the number-conserving Hamiltonian 

N 

77 = - ^ - tb^bj+i IF6]4+i%%+i + h-c. 

i=i 

(18) 

In [T7], evidence of a topological phase in this system 
was found, in the large N limit, at around the 1/3 fill¬ 
ing fraction, and near ta + tb = W. In this region, 
the two ground states in the even and odd a-wire par¬ 
ity sectors have an energy gap that grows exponentially 
small in the system size m- In figure [^top) , we plot 
exact results of the number difference 0 i in this phase 
(ta = tb = 0.5IT, n « 1/3), as a function of system size 
{2N). This shows signs of exponential decay, however 
finite size effects cause non-linearities in the curve. 

As the particle number is the same in both of the states 
we are interested in, our probe cannot be al. Instead, an 
excitation must involve hopping particles between the a 
and b wires. Thus, we propose pij = aJ4 as a potential 
probe of the system. If our Majorana bound states are 
localised at the edges of the wire, we would expect pij = 
|(^£;|Py|'I'o)| to be approximately 0 unless both i and 
j are near the edges of the system. Unfortunately, the 
small system sizes we use here obscure these signatures. 
In figure [^bottom), we plot J2j Pij and PhnE where 
hff is taken from the middle of the system with system 
size N (other parameters are ta = h = 0.5IU, n « 1/3). 
We see that the probe density in the middle of the system 
is much smaller than at the edges. However we do not 
have enough data to be sure that one survives whilst 
the other drops to zero in the large N limit. As DMRG 
studies have been performed on this system on up to 
A = 36 [T7], definite results should be obtainable. From 
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the above discussion, we expect the Majorana operators 
to take the form ^ a\bjN-\-h.c. with i and j on the same 
ends of the wire, and N a product of number operators. 

The results derived in this paper do not present a com¬ 
pleted study of an MBS in the many-particle framework. 
It is possible to immediately include larger numbers of de¬ 
generacies; with 2" ground states we will obtain 2n MBSs 
by repeating our procedure on pairs of ground states. 
However, other important questions remain unanswered, 
like the exact nature of braiding, fusion and measure¬ 
ment in this formalism. We hope these concepts could 
be inherited naturally from existing work. 
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Proof that localisation condition is exact 


Our proof proceeds by giving a method to construct 
the Majorana operator when the localisation condition 
is satisfied, and proving that the method works consis¬ 
tently. Let us assume that we have satisfied the localisa¬ 
tion condition for J. Then, for all J C we can define 
nj ■■= {'SE\nj\'i' e) = ('I'oI^jI^'o)- Then, let A = O or 
E, and define 


:= 


(nj - nj) I^'a). 


y/ nj — vtj 

Then, fix AT, J C J, and we can calculate 


(19) 


flK\^A,j) = 


- 2 {"njuK — njfiK) I'I'a) 

nj - Uj 


1 


{nKuJ - ,KUj) 


flj - Uj 

-njiriK - n\Y/‘^\^A,K) + (riifuj - nKnj)\^A) 


Thus, the action of all fiK operators remains in the space 
S spanned by {['I'a), for J S J' and A = 0,E. 

Then, we can calculate the matrix elements for I, J,K C 

J 


{^I,A\nj\^K,A) 


{niujuK — ninjuK — nnjjn-K + ninjUK) 
^J(ni - nj){nK - n\) 


Importantly, these are independent of A, and so within 
the subspace 5, our number operators act as the block di¬ 
agonal matrices fij = nj®\. Although the |5'yA) states 
are not orthonormal, this fact will not change under the 
Gram-Schmidt procedure applied equally to both sectors 
(as the overlap {'^.j.a\^k,a) is again independent of A). 

The number operators hj (note small j) are mutu¬ 
ally commuting, and so they can be simultaneously di- 
agonlised within this subspace S. Within this eigenba- 
sis, we can assume that all eigenvectors are separated by 
eigenvalues of at least one number operator, except that 
they must come in pairs of opposite parity (due to the 
hj = fij ® structure of the number operators). Each 
of these pairs of eigenvectors define a two-dimensional 
subspace Hi. 

The most general form our Majorana operator can 
take (to commute with each hj) is then 7 = ^ 

with erf = cos(()))crf -I- sin((/))crf. The coefficients 7 ^ are 
fixed by transforming back to the original basis (of gram- 
schmidt orthonormalised {['I'a), I'h^.j)}), and insisting 
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that 7 contains no terms of the form |'I'y 4 )(‘hB.j|- In 
matrix form, the transformation looks as follows: 



( 20 ) 


Here, N is the dimension of S divided by 2 (which has a 
maximum of 2 l‘^l, but does not necessarily achieve this). 
Importantly, we see that there are N degrees of freedom 
in our Majorana, and N conditions fixed by the top row 
of the original basis. This counting shows that as long as 
the conditions are not contradictory, at least one solution 
will exist. 

Let [/ 0 1 be the unitary matrix that transforms from 
the fij eigenbasis to the original basis. Then, our condi¬ 
tions can be written in the form 

= Sik- ( 21 ) 

3 

If we write 7 = ( 71 ... Jn), and define the matrix Aij = 
U*iUjk, than we have reduced our problem to solving the 
set of linear equations H 7 = (1,0,..., 0). A solution to 
this will exist as long as the first row of A is not a lin¬ 
ear combination of the others. To see that this is always 
the case, note that the row vectors in A are obtained 
by multiplying the orthogonal column vectors in [/ by a 
diagonal matrix Dj^k = This is a shearing op¬ 

eration, which cannot transform orthogonal vectors into 
linear combinations of each other unless it projects them 
to the zero vector. And, the first row of A is not pro¬ 
jected to the zero vector as it is a termwise square of the 
elements of the first column vector in U, which has norm 
1. Thus, it is always possible to construct at least one 


localised Majorana operator via the above method when 
the localisation condition holds. 

From the above proof, it is possible to calculate the 
localised Majorana operator when the localisation condi¬ 
tion is explicitly satisfied, by construction of the matrix 
A. In the case of small numbers of sites, we can obtain a 
closed form of the coefficients of the Majorana operator. 
Let us write 


1 I ’ 

1 - ni - n2 -I- ni2 
^ ^ (n,-l)(ni-nf) 

(ni2 - ni)(l - ni - 712-I-ni2) ’ 

^ ^ (tT, - l)A/(7r, -Tlf) (7712 - 77 ^ 2 ) 

( 7^12 - 77j)(l - 77i - 772 -I- 7712 ) ’ 

^ _ 77i772(1 -771 - 772 -1- 7772 ) -b 77i2(77i772 - 77i2) 

(7712 - 772 ) (1 - 77i - 772 -b 7712 ) 

Here, 7 = 1 or 2, and i is the opposite site to i. Note that 
this solution assumes the linear independence of all the 
I'I’a.j); when this is not the case, some of these denom¬ 
inators will be infinite, and a different set of solutions 
will need to be calculated. These coefficients were used 
in the construction of the Majorana operators for the 
double quantum dot. 


7 = 7b'+ ^ {e^’^Gj,K\^o,j){^E,K\+e ^‘^Gj^k\^e,k){^o,j\) , 


and then solving the linear equations above will give us 
the coefficients of the Gjk- To remove a single site, 
we have only one coefficient to solve for, and we obtain 
Gip = 1. The solutions for two sites are as follows (with 
Gj,k = Gk,j) 


Ai o = 


''(771 - nl){n2 - nl) 



